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ABSTRACT 

A periodic variation in the pulse timings of the pulsating hot subdwarf B star CS 1246 was recently 
discovered via the O-C diagram and suggests the presence of a binary companion with an orbital 
period of two weeks. Fits to this phase variation, when interpreted as orbital reflex motion, imply CS 
1246 orbits a barycenter 11 light-seconds away with a velocity of 16.6 km s _1 . Using the Goodman 
spectrograph on the SOAR telescope, we decided to confirm this hypothesis by obtaining radial 
velocity measurements of the system over several months. Our spectra reveal a velocity variation with 
amplitude, period, and phase in accordance with the O-C diagram predictions. This corroboration 
demonstrates that the rapid pulsations of hot subdwarf B stars can be adequate clocks for the discovery 
of binary companions via the pulse timing method. 

Subject headings: stars: oscillations - stars: binaries - stars: subdwarfs - techniques: radial velocities 
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1. DETECTING BINARIES WITH PULSE TIMINGS 

Binary star systems are significant in astronomy be- 
cause their orbits often permit the calculations of im- 
portant astrophysical parameters including mass, radius, 
and density. Binarity plays an apparently crucial role in 
the story of the hot subdwarf B (sdB) stars. Modeled 
as extended horizontal branch stars in globular clusters, 
sdBs are presumed to have He-burning cores surrounded 
by a thin layer of Hydrogen. Models show they are the 
progeny of red giant branch stars that somehow lost the 
majority of their outer Hydrogen layer before reaching 
the Helium flash. Many theories attempt to explain the 
expulsion of this mass, and most of them invo ke the pres- 
ence of a binary compa nion, either stellar ( Ha n et al.l 
120021 12003D or substellar (|Sokerlll998f) . Observations sup- 
port this picture; several studies h ave revealed a large 
binary fraction amongst the sdBs ([Maxted et al.l 120011 : 
ICopperwheat et all 1201 lh . Constraining the orbital pa- 
rameters of these systems may shed light on the evolu- 
tionary histories of hot subdwarf stars, and many studies 
are being carried out w ith this goal in mind (jSchuh et al.l 
120101: iGeier et al.ll2011l ). For further detail s on hot sub- 
dwarfs, we refer the reader to iHeberl (|2009[) . 

The majority of hot subwarf orbital parameters cal- 
culated to date come from spectroscop ic radial velocity 
(RV) measurements (|Geier et al.l 1201 lL Table A.l). In 
most sdB binary systems, the companion contributes a 
negligible amount of the system light, so measurements 
are confined to reflex motion of the sdB. Although the 
RV method is quite effective at uncovering reflex mo- 
tions, the technique is heavily biased towards systems 
with small separation distances and high masses. 

For binary systems with wider orbits, the light-travel 
time (LTT) techniqu^3 provides a viable alternative to 
the RV method. As each member in the system orbits the 
barycenter, the line-of-sight distance from it to the Earth 
oscillates; thus, light from the object sometimes arrives 
later than usual and other times earlier. If the binary 



contains a pulsating star, this changing light-travel time 
manifests itself as a phase oscillation in the observed mi- 
nus calcul ated (O-C) diagram constructed from the pulse 
timings. ISilvotti et al.l (|2007D used this technique and 
found evidence of a 3.2-Mj planet orbiting the sdB star 
V391 Peg, which demonstrated that some planets can 
survive the red giant expansion of their host stars. Sev- 
eral reviews have been published describing the proce- 
dures, shortco mings, and successes of this method (e.g., 
ISterkenl[2fJ05l) . 

Few systems exist with orbital parameters easily de- 
tectable by both the RV and pulse timing methods. In 
fact, out of all types of pulsating stars, only two have 
been observed with RV measurements in accordance with 
the result s of the LTT technique: AW Persei, a cl assical 
Cenh eid (IVinkd Il99l and S Z Lyn, a 6 Scuti jjSzeidl 
[19831: IBardin fc Imberti Il98i IMoffett et all 11988ft . No 
such corroboration has been demonstrated yet using pul- 
sations with periods on the or der of minutes, as exist 
in the pulsating w hite dwarf ([Montgomery! l2009f ) and 
hot subdwarf stars (|0stensenl 120091 ) . Nonetheless, much 
time and effort is being spent accumulating phase mea- 
surements of these types of pulsators in hopes of dis- 
covering binary members down to pl anetary-size (e.g., 
ISchuh et al.l[2010t iMullailv et al.ll2008l) . 

Recently, we discovered fortnightly phase oscillations 
in the pulse ti mings of the sdBV r star CS 1246 via the 
O-C diagram (B arlow et al.l l2Ol~0) . An orbiting body 
with minimum mass of 0.13 Mq provided the simplest 
explanation for the phase wobble, but since no secondary 
light could be detected from the system, this hypothesis 
remained unconfirmed. Here, we present corroborating 
radial velocity evidence for a binary companion around 
CS 1246. Our result confirms that the rapid pulsations 
of hot subdwarf B stars can be sufficient clocks for the 
detection of companions using the O-C diagram. 
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1 also known as the pulse timing method when pulsations are 
used as the clock 
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Fig. 1. — SOAR/Goodman spectra of CS 1246. Top panel: Example spectrum used to determine the RV of CS 1246 on a single night. 
The spectrum shown, which is a combination of five 371. 7-s exposures, has a resolution of 0.52 A. Bottom panel: Average of all spectra 
used in this study. The spectra were individually velocity-shifted to remove the RV variation from the orbital reflex motion. The average 
signal-to-noise ratio is approximately 270 per resolution element. 



2. THE PULSE TIMINGS OF CS 1246 

Using the 0.41-m PROMPT array, we aggressively 
monitored the pulsations of CS 1246 for four consecu- 
tive months in 2010 to look for secular variations in the 
pulse timings. The phase of the 371. 7-s pulsation mode 
(the only confirmed mode) was used to construct an O-C 
diagram, which revealed a parabolic secular trend with 
overlying sine wave. For full details on these observa- 
tions, including the construction and inte rpretations of 
the O -C diagram, we refer the reader to iBarlow et al.l 
(l2T?Tl . 

A secular decrease in the pulsational period on the or- 
der of 1 ms every 1.7 years can explain the parabolic 
component. We must note, however, it is plausible that 
the quadratic term is instead a sine wave with period 
much longer than our baseline, which could masquerade 
as a parabola until a wider phase coverage is obtained. 

The sinusoidal trend is undoubtedly more convincing. 
Having undergone 30 cycles over the course of our ob- 
serving baseline, the phase oscillation appeared in both 
the amplitude spectrum of the combined light curve (as 
a splitting of the main mode) and in a Fourier transform 
of the O-C diagram itself. Several phenomena can be 
invoked to account for this variation, but the presence 
of a previously-unknown binary companion provided the 
simplest explanation. As the orbiting body tugs on CS 
1246, the light-travel distance changes along our line-of- 
sight, resulting in a phase oscillation of the pulsations. 
This hypothesis is easily falsifiable since it makes the fol- 
lowing prediction: if an orbital reflex motion is the cause 
of the phase oscillation, CS 1246 should exhibit line-of- 
sight radial velocity variations with a semi- amplitude of 
16.6 km s _1 and period of two weeks. Moreover, the 
variation should be 90 degrees out-of-phase with the os- 
cillation in the O-C diagram. Table 1 presents the full set 
of orbital para meters predict e d from the pulse timings, 
as published bv IBarlow et al.l (|2011( ). 

3. SPECTROSCOPY 

We began taking spectra of CS 1246 in January 2011 
using the Goodman spectrograph on the 4.1-m SOAR 
telescope. Appendix A presents the details of our ob- 
servations. Since an RV oscillation with an observed 



TABLE 1 

Orbital Parameters from the O-C Diagram 



Par am 


Value 


Error 


Units 


Comments 


n 


14.105 


± 0.011 


days 


orbital period 


Ki 


16.6 


± 0.6 


km s _1 


CS1246 RV semi-amp 


f 


0.0066 


± 0.0007 


M© 


mass function 


a 


0.0963 


± 0.0003 


AU 


separation distance 


msini a 


0.129 


± 0.005 


M 


companion mass 



a assuming the canonical sdB mass 0.47 Mq for CS 1246 

semi-amplitude of 8.8 km s " 1 is generated by the pul- 
sational motions themselves (jBarlow et al.l[2010h . we set 
the exposure time of each spectrum to the pulsation pe- 
riod, 371. 7-s. This technique washes out the pulsational 
RV variations and makes it easier to detect changes in 
the systemic velocity. We used a 0.46" slit and a 2100 
mm VPH grating from Syzygy Optics, LLC, to cover 
the spectral range 3690-4390 A with a resolution of 0.52 
A (0.17 A per pixel). To avoid contaminating the spec- 
tra with light from a dim star a few arcseconds to the 
West, we set the position angle along the slit to 270 de- 
grees East of North and sacrificed observing at the par- 
allactic angle. Flexure in the spectrograph as an object 
is tracked changes the wavelength solution slightly over 
several hours; consequently, we avoided taking more than 
30 minutes of stellar spectra in a row without obtaining 
comparison spectra of an Fc lamp. When time permit- 
ted, we als o observed a n A-type radial velocity standard 
(HR 3383; iFekell I1999D to determine the velocity zero- 
point. In total, we took 149 individual spectra of CS 
1246 over a six- month period. 

Each spectral frame was bias-subtracted and flat- 
fielded in IRAF. Additionally, we subtracted "pseudo" 
darks to remove a coherent source of stray light in the 
spectrograph with a count rate of nearly six photons per 
minute per pixel. These calibration frames were obtained 
by placing the spectrograph in the same configuration 
used for the target observations and integrating with the 
shutter open (in darkness, with the flip-mirror of the tele- 
scope closed). We did not flux-calibrate the spectra as 
time did not permit us to observe spectroscopic standard 
stars each night. 

We combined together all spectra taken consecutively 
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on a given night to produce "master" spectra. On most 
nights, we were able to produce two master spectra, with 
total integration times around 26 minutes and signal-to- 
noise ratios of 60 per resolution element, on average. Fig- 
ure [T] presents an example of one of these spectra plotted 
above the mean spectrum created using all of our data. 

4. THE RADIAL VELOCITY CURVE 

We employed two methods to measure the radial ve- 
locities in our master spectra. Both are iterative pro- 
cesses during which we (i) create a mean spectrum from 
all individual spectra, (ii) measure velocities of individ- 
ual spectra based off of information from the mean, (iii) 
create a new mean spectrum by velocity-shifting the in- 
dividual spectra, and (iv) repeat the above steps until 
the solution converges. Further details on these methods 
follow. 

(1) Lorentzian+Gaussian Fits to H Ba lmer Lines: 

Using the MPFIT routine in IDL IMarkwardtl 
(2009), which employs the Levenberg-Marquardt 
method, we first fit the sum of a Gaussian and 
Lorentizian to the individual H Balmer line profiles 
in the mean spectrum. Fixing the shapes of the 
profiles to those fit in the mean spectrum, we 
then fit the Balmer line profiles in the individual 
spectra and used the best-fit centroid values 
to compute velocities. The mean velocity of 
each master spectrum was calculated from the 
weighted-average of the Balmer line velocities. 
We created an improved template by de-shifting 
the individual spectra using the previous RV 
results and repeated the process until the solution 
converged. Balmer lines beyond H9 were ignored 
since the low signal-to-noise ratio in this part of 
the spectrum resulted in inconsistent fits. We esti- 
mated the velocity error in each master spectrum 
using the standard deviation of the individual 
line-profile measurements. If a night contained 
more than one master spectrum, we averaged the 
resulting measurements. To check the zero-point, 
we also measured the radial velocity of the velocity 
standard HR 3383. As our measurement agreed 
with the published value to within the errors, 
we did not apply any offsets to our measured 
velocities for CS 1246. 

(2) Cr oss-Correla t ion: Following a procedure similar 

to lSaffer et al.l ([20011 ). we used the fxcor routine in 
IRAF , which uses the technique of iTonrv fc Davisl 
([1979D . to compute the velocity shift of each master 
spectrum. All spectra were combined into a tem- 
plate that was cross-correlated against the individ- 
ual, unshifted spectra. An improved template was 
then created by de-shifting the individual spectra 
using the previous cross-correlation results, and the 
process was repeated until the solution converged. 
We converted the relative velocities to absolute ve- 
locities by cross correlating the template with the 
spectrum of our observed radial velocity standard, 
HR 3383. Errors were taken from the output of the 
fxcor task in IRAF. If a night contained more than 
one master spectrum, we averaged the resulting RV 
measurements. 



After determining the velocities, we used the IRAF task 
rvcorrect to calculate the heliocentric corrections at our 
observing epochs and applied these corrections to our 
measurements. Results are shown in Appendix A. 

RV curves produced using both methods display sim- 
ilar variations. Figure [2^ presents the curve from 
the Lorentzian+Gaussian method as an example. The 
dashed line marks the variation predicted by the O-C di- 
agram, with the overall systemic velocity fit as a free pa- 
rameter. Figure[2L shows the same measurements phase- 
folded on the 14.1-day periodicity observed in the O-C 
diagram. 

Although inspection of these figures makes clear the 
agreement between data and predictions, we performed 
non-linear, least-squares fits of sine waves to the RV 
curves to quantify the similarities, leaving all param- 
eters free. Table 2 presents the best-fit parameters 
above the O-C diagram predictions; the dotted lines 
in Figure [5] show the best-fitting sine wave to the 
Lorentzian+Gaussian fit RV points. The period and am- 
plitude calculated using both methods agree with those 
predicted by the O-C diagram. Our RV data allow us to 
calculate for the first time the system velocity. By fitting 
a sine wave to the data with the period, amplitude, and 
phase fixed to their O-C calculated values, we derive a 
systemic velocity of 67.3 ± 1.7 km s _1 . 

5. CONCLUSIONS 

We have presented radial velocity data confirming the 
binary nature of CS 1246. The period, amplitude, and 
phase of the RV variation agree with those predicted by 
the O-C diagram when interpreting the observed phase 
oscillation as an orbital reflex motion effect. Although 
this is not the first instance such an agreement between 
pulse timing and velocity results has been noted, this 
study is the first using pulsations on the order of min- 
utes and the first for a hot subdwarf star. Our result 
demonstrates that the rapid pulsations of sdB stars can 
be adequate clocks for detecting binary companions us- 
ing the light-travel time technique. 

Since the error bars of the orbital parameters calcu- 
lated from the O-C diagram fits are smaller than those 
calculated from the RV curve, our characterization of the 
companion remains the same as originally presented in 
iBarlow et all (|2011[ ) (see Table 1). Using the canoni- 
cal sdB mass (0.47 M ) for CS 1246, we calculated the 
minimum mass of the companion to be 0.13 M©. Since 
no eclipses or other photometric indications of the com- 
panion have been observed, the orbital inclination angle 
remains unknown, and we are limited to computing its 
minimum mass. There is a 95% probability the com- 
panion mass is less than 0.45 M©, assuming an isotropic 
distribution of inclination angles. Thus, CS 1246 is most 
likely orbited by a low-mass white dwarf or M-type main 
sequence star. 

According to the bina ry population synthesis models of 
lHan et all ([20021 . [20031 , an sdB binary with a two- week 
orbital period could be the product of both common en- 
velope ejection and stable Roche lobe overflow formation 
channels. The latter is unlikely for the CS 1246 system 
since it predicts an early-type main sequence companion, 
signatures of which would be visible in our spectra but 
are not. Both the "first" and "second" common enve lope 
ejection channels presented in lHan et al.l (|2002t [20031) are 
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Fig. 2. — Heliocentric radial velocities of CS 1246. (a) Top panel: RV measurements derived from Lorentzian+Gaussian fits to the H 
Balmer absorption-line profiles. The dashed line marks the velocity curve inferred from the phase oscillation in the O-C diagram, under 
the assumption it is caused by orbital reflex motion. Note the agreement in the period, amplitude, and phase. The dotted line shows the 
best-fitting sine wave to the data, with all parameters left free. Bottom panel: Residuals after subtracting from the data the RV curve 
predicted by the O-C diagram, (b) Top panel: RV curve folded on the period predicted by the O-C diagram (14.105 days), plotted twice 
for better visualization. The dashed and dotted lines show the O-C diagram prediction and the best-fit sine wave, respectively. Bottom 
panel: Residuals after subtracting from the data the RV curve predicted by the O-C diagram. The standard deviation is 4.7 km s _1 



more plausible. In these scenarios, the sdB progenitor, a 
red giant star, transfers mass dynamically to its compan- 
ion prior to reaching the tip of the giant branch. Even- 
tually a common envelope forms, the orbit shrinks due 
to friction and tidal interactions, and the envelope of 
the giant star is ejected, leaving behind a close binary 
composed of an sdB star with a main sequence ( "first" 
channel) or white dwarf ("second" channel) companion. 
Their models predict similar sdB masses for both com- 
panion types. If the companion is an M-dwarf, the sys- 
tem might show an infrared excess, but looking for this 
photometric signature is complicated by the presence of 
the Coalsack Dark Nebula, which reddens the light sig- 
nificantly. Due to these complications, we are left with 
no choice but to continue improving measurements of the 
orbital parameters of this system using additional pulse 
timing data. 

We acknowledge the support of the National Science 
Foundation, under award AST-0707381, and are grateful 
to the Abraham Goodman family for providing the finan- 



TABLE 2 

Velocity Curve Fit Parameters 



Method 


Period 


Amplitude 


Systemic RV 




(days) 


(km s" 1 ) 


(km s" 1 ) 


Lorentz+Gauss 


14.03 ± 0.06 


19.8 ± 2.7 


67.5 ± 1.9 


Cross-Corr. 


13.92 ± 0.15 


13.9 ± 2.9 


66.9 ± 2.7 


O-C Diagram 


14.104 ± 0.011 


16.6 ± 0.6 


no prediction 
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APPENDIX 



TABLE Al 

Log of Master Spectra & Velocities 



Date" 


HJD a 


Integration Time 


Radial Velocity 6 


(UT) 


(HJD) 


(<0 


(km 


s 




2010-11-09 


2455509.85674 


4 x 371.7 


76.5 


± 


4.9 


2011-01-21 


2455582.85745 


7 x 371.7 


49.7 


± 


4.2 


2011-01-22 


2455583.83098 


8 x 371.7 


51.1 


± 


2.6 


2011-01-23 


2455584.83992 


10 x 371.7 


50.7 


± 


4.0 


2011-01-24 


2455585.83775 


10 x 371.7 


53.7 


± 


4.2 


2011-01-25 


2455586.84322 


9 x 371.7 


49.1 


± 


5.0 


2011-01-26 


2455587.76537 


16 x 371.7 


61.8 


± 


3.4 


2011-01-28 


2455589.83534 


11 x 371.7 


76.9 


± 


3.2 


2011-02-07 


2455599.78876 


6 x 371.7 


54.8 


± 


2.2 


2011-02-16 


2455608.85764 


6 x 371.7 


71.6 


± 


2.6 


2011-02-17 


2455609.79579 


8 x 371.7 


63.8 


± 


3.2 


2011-02-18 


2455610.86934 


8 x 371.7 


58.5 


± 


2.7 


2011-02-19 


2455611.80540 


12 x 371.7 


46.2 


± 


1.4 


2011-02-23 


2455615.86702 


5 x 371.7 


69.6 


± 


3.5 


2011-02-24 


2455616.84545 


6 x 371.7 


79.7 


± 


4.6 


2011-03-18 


2455638.51437 


10 x 371.7 


55.3 


± 


5.1 


2011-03-22 


2455642.63103 


7 x 371.7 


48.8 


± 


2.9 


2011-05-24 


2455706.48698 


3 x 371.7 


81.9 


± 


7.8 


2011-05-26 


2455707.72320 


3 x 371.7 


63.0 


± 


4.0 



° at mid-integration 

b heliocentric velocities from Lorentzian+Gaussian fitting method 



